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The critical behavior of a model describing phase transitions in 3D antiferromag- 
nets with 2A^-component real order parameters is studied within the renormalization- 
group (RG) approach. The RG functions are calculated in the three-loop order and 
resummed by the generalized Pade-Borel procedure preserving the specific symmetry 
properties of the model. Anisotropic stable fixed point is found to exist in the RG 
flow diagram for N > 2 and lie near the Bose fixed point; corresponding critical 
exponents are close to those of the XY-model. The accuracy of the results obtained 
is discussed and estimated. 
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In this Brief Report, we study within the field-theoretical RG approach in three dimen- 
sions (3D) the critical behavior of a model with three quartic coupling constants describing 
certain antiferromagnetic and structural phase transitions. The Hamiltonian of the model 
reads 



H = j d X -{m^ipaiPa + V^cN^a) + -^Wa^a) + -^^^a + -^^2(5-1^2(5 



(0.1) 



where ipa is a real vector order parameter field, a = 1, 2, ... , 2N , [3 = 1,2, . . . , A^. For 
Uq = Hamiltonian ( p.l|) describes non-interacting anisotropic XF-models, for Zq = 
it reduces to that of the cubic model. When N = 2 the expression ( |0.1|) is relevant to 
the structural phase transition in the Nb02 crystal and, for vq = Zq, to antiferromagnetic 
transitions in TbAu2 and DyC2- Another physically interesting case A^ = 3 corresponds to 
antiferromagnetic phase transition in the K2lrClQ crystal and, for vq = zq, to those in TbD2 
and Nd [§. 

This model was studied earlier by the e-expansion technique and directly in 3D within 
the lower (one- and two-loop) perturbative orders 0]^. The main result of those investiga- 
tions was the existence of the nontrivial (anisotropic) stable fixed point for N > 2 which 
gives rise to a new universality class with a certain set of critical exponents. The lower-order 
approximations, however, are known to lead to crude quantitative and, sometimes, to con- 
tradictory qualitative results, especially for systems with complicated symmetry B . That is 



why recently the critical behavior of the model ( p.l|) has been analyzed in the third order in e 
1^. Investigation of the fixed points stability and the calculation of marginal dimensionality 
Nc of the order parameter separating two different regimes of critical behavior confirmed 
that the model ( |U.1| ) possesses the anisotropic stable fixed point for physically interesting 
cases N = 2 and = 3. However, the eigenvalue exponents for this point turned out to be 
twofold degenerate in the one- loop approximation . Such a degeneracy decreases markedly 
the accuracy expected within a given approximation and makes the calculating the eigen- 
value exponents in higher orders in e very difficult, while resummation of shorter series fails 
to provide proper numerical estimates. It is reasonable therefore to study the model ( p.l|) 
by means of some alternative approach. Below, the critical behavior of this model will be 
investigated by the RG technique in 3D in the three-loop approximation, aiming to improve 
existing qualitative and numerical results. 

The character of the critical behavior is determined by the RG equations for quartic 
coupling constants. Calculating the /3-functions entering their right-hand sides and critical 
exponents 7 and rj as functions of the dimensionless coupling constants u, v, and z in the 
three-loop approximation, we obtain: 
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I3u = u\l - u- jrriiSv + z) + 



8(Ar+4)3 



2(41 + 95)^2 + SOOuv + lOOuz 

3 



27(Af+4)2 

"(5.39577 + 109.88075Ar + 199.64042)^^ 
+ (39.88127 N + mM155)u^v + (13.29376 N + 164.61385^^ 

+ (3.73134 N + 302.86778)Mf 2 + 58.35119m;^ + (1.24378 N + 81.50553) 



X uz 



2 + 65.93728^3 + 3.73134^^^ + 21.97909^^^ ^ 9.35389^3" 



(3v ^ V- j^(6uv + + ^z^) + ^f^^^ [2(23 N + 185)u''v + 62Auv^ 



-5.00443 



^Z^ 



+ A6uvz + 5Auz'^ + 231v^ + 23vz'^ + 18z^] - ^^j^^ 
+ 83.70780 N + 469.33398)^3^; + (4.47812 N + 1228.60593)^2^; 
+ (-7.50664 N + 115.06965)ii2-i;z + (2.99978 N + 98.15522)w 
+ 957.78168mv3 + 4:A7812uv^z + 135.72564mv^2 ^ 60.68074^^3 
+ 255.92974^^ + 4.10473vV + 49.83176^;;s3 + 5.05072;s4], 

(^^ = -^{l - jvT4(6^* + 3^ + 2^) + 27iN+4r [2(23 N + 185)^2 + 462m^; 

-5.00443 + 83.70780 



+ 262uz + 69v^ + 162vz + Alz^^ sUv+W 
+ 469.33398)^3 + (-4.52123 + 934.14027)^2^; + (4.49249 
+ 507.69053)ii2^ + 411.65505^^2 + 550.60474iivz + 196.40638^^2 



+ 65.93728t;3 + 108.78727t;2^ + 118.33643t;^2 ^ 21.85596^3 



7-^ = 1 
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4(Ar + 4) 
0.04813 



2(A^+ 1)m + 3^; + ^ 



1 



(0.2) 



2(Ar + 1)k^ + 6uv + 2uz + ?>v^ + z^ 



4(Ar + 4)2 

4(Ar + 1)2^3 + isi^jv + l)M2-t; + 6(A^ + l)u^z + 6(A^ + A)uv'^ + I2uvz 

0.06182 



(A^ + 4)3 

+ 2(Ar + 2)uz'^ + 9^3 + 3^2;^ + 2,vz'^ + z^ 



4(A^ + l)(Ar + 4)M3 



(A^ + 4)3 

+ 18(Ar + 4)^2-1; + 6(Ar + 4)^2^ + 81hv2 + ISuvz + 21k^2 ^ 27^3 + Qvz^ + 4^31 ^ (o.3) 



27(Ar + 4)^ 



2{N + 1)^2 + Quv + 2uz + 3^2 + z^ 



+ 



0.00309 
(AT + 4)3 



4(Ar+ l)(iV + 4)M3 



+ 18(A^ + A)u^v + 6(A^ + 4)^2^ + 81mv2 + IMvz + 21m^2 ^ 27^3 + Qvz"^ + 4^3 _ (o.4) 
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To extract the physical information from these divergent series, the Pade-Borel resumma- 
tion procedure will be used. Since the expansions of quantities depending on three variables 
u, V, and z are dealt with, the Borel transformation is taken in a generalized form: 

fiu, v,z) = Y: c,,,uVz'' = I e-'F{ut, vt, zt)dt, F{x, y,w) = Y: (-'f-^l^r ^^'^^ 

ijk Q ijk y''' ' 3 ' )■ 

To perform an analytical continuation, we address to the resolvent series 

F{x, 2/, ^, A) = 5: A" E E ''"^'""'""^ f (0.6) 

n=0 1=0 m=0 ^• 

which is a series in powers of A with coefficients being uniform polynomials in u, v, z, and 
then use Fade approximants [L/M] in A at A = 1. The approximant [3/1] is employed 
for resummation of three-loop RG series Eqs. ( p.2| ). The coordinates of all fixed points 
are determined for two most interesting cases N = 2 and = 3. The results obtained 
are presented in Table I, which contains also, for comparison, analogous estimates found 
earlier from the two-loop RG series The global structure of the RG fiows in the three- 
loop approximation is shown in Fig.l (see the printed paper: Phys. Rev. B 59 (1999) 
8363-8366). 

The three-loop contributions to the /3-functions are seen to affect the location of the 

fixed points considerably. Of special interest is the fortune of the nontrivial, anisotropic 

fixed point 8. In the one-loop approximation for = 2 it coincides with the Heisenberg 

fixed point 2 In higher orders this degeneracy is lifted out and the point 8, remaining 
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a stable one, moves toward the Bose fixed point 5 lying within the plane v = z. With 
increasing A^, the former becomes closer and closer to the latter, indicating the tendency of 
the 0(2A^)-symmetric interaction to vanish at criticality. This resembles the behavior of the 
cubic model that is known to split into n non-interacting Ising models when n oo. 

The resummation procedure changes the perturbative expansions for /^-functions by the 
complicated non-polynomial expressions. Does this procedure preserve the symmetry prop- 
erties of the system? The model ( p.l|) possesses, apart from usual, the special symmetry 
properties which have been already used for testing a validity of approximations employed . 
For example, if the field ipa undergoes the transformation v?2/3-^i (v'2/3-1 + ^213), f2i3 

{'^2^-1 — V213), the coupling constants are also transformed: 

1 1 
u^u, v-^-{v + z), z^-{3v-z) , (0.7) 

but the structure of the Hamiltonian itself remains the same. Since such a transformation 
does not affect the RG equations, it can, at most, rearrange the fixed points leaving their 
coordinates unchanged PJ^,§]. To check up whether such an invariance holds for the fixed 
points found, let us apply the transformation ( p.7| ) to the content of Table I. Then points 
1, 2, 5 and 8 will stay at their places while the coordinates of points 3 and 4 will turn into 
those of points 6 and 7 respectively and vice versa with accuracy of order of 10^"^. Since the 
fixed point coordinates were evaluated numerically with just the same accuracy, it means 
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that our resummation procedure exactly reproduces the symmetry properties discussed. 

Let us determine further the critical exponents. The coefficients of the series for rj 
rapidly diminish, therefore its value may be found by direct substitution of the fixed point 
coordinates into the expansion (|U.4|). In the case of susceptibility exponent corresponding 
series ( p.3|) have to be resummed. Using the generalized Pade-Borel procedure, we evaluate 
the exponent 7 and then, addressing the values of rj, calculate critical exponent a by the 
scaling relations. The results for N = 2 and = 3 are presented in Table II. As is seen, 
the critical exponents of the fixed point 8 only slightly differ from those of the Bose one 
reflecting a closeness of both fixed points in the (u, v, z) space. 

The critical exponents for the fixed points 2 (Heisenberg), 5 (Bose), and 3 (Ising) thus 
found are close to their high-precision analogs resulting from the six-loop RG expansions 
P]; the differences do not exceed 0.02-0.03. For the anisotropic fixed point 8 the situation 
is more complicated. In the three-loop approximation it is three-dimensionally stable and 
characterized by negative value of the specific heat exponent a, while the Bose fixed point 
looks unstable since it has a positive a. In fact, however, the XF-like critical behavior in 3D 
is described by the negative a as is known both from the six-loop RG calculations @] and 
recent extremely accurate measurements |p . Moreover, a sign of the exponent a at the Bose 
fixed point determines either the 0(2A^)-symmetric interaction is relevant near this point or 
not ||lO|: if a < 0, the Bose fixed point should be stable with respect to this interaction. It 



means that what is really stable is the fixed point 5 while the stability of the fixed point 8 
within the plane v = z is an artifact of the three-loop approximation. 

Keeping this in mind we can estimate an accuracy of the anisotropic fixed point coor- 
dinates given by the three-loop RG series. Because of the obvious topological reasons (see 
Fig.l), the point 8, being in fact unstable, should have Uc < 0. At the same time, since the 
true value of a at the Bose fixed point is very small p],^ the points 5 and 8 should be very 
close one to another and the modulus of Uc for the point 8 should be very small as well. 
Hence, as can be deduced from the numbers in the last column of Table I, the three-loop 
approximation predicts locations of this point for N = 2 and = 3 with errors about 0.2 
and 0.1 respectively. 

Dealing with the theory without a small parameter and the short (three-loop) perturba- 
tive series one would refer to such an accuracy as satisfactory. On the other hand, numerically 
small errors lead in this case to qualitatively incorrect results making the situation rather 
unfavorable. The point is that in 3D the model ( |0.1| ) is almost identical to some marginal 
system for which a = at the Xy-like criticality and, as a consequence, fixed points 5 and 
8 coincide. Since such a marginality is far from to be seen in the one-loop approximation, 
it hardly manifests itself within a perturbation theory, even in higher orders. 

The " near- marginality" is not a unique feature of the model ( p.lD being typical for 3D 

systems with several coupling constants. Thus, for the cubic model the marginal value of n 
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that separates two different regimes of critical behavior is numerically close to the physical 



value n = 3 [0. The model describing phase transitions into chiral states has the marginal 
dimensionality of the complex order parameter Nc2 (separating domains of continuous and 
first-order transitions) that is also very close to the physical value = 2 0|, etc. Hence, to 
answer the question about the type of the critical behavior of such "unconvenient" models 
the higher-order RG analysis should be carried out. 

The calculation of the RG functions of the model ( |0.1|) in the next, four-loop approxi- 
mation is an extremely difficult problem. The point discussed, however, may be cleared up 
without performing additional RG calculations. Indeed, let us trace how the numerical value 
of the exponent a for the Bose fixed point depends on the order of the RG approximation. 
The estimates for a obtained in the one-, two-, three-, four-, five-, and six-loop orders using 
the Pade-Borel resummation technique are 0.125, -0.012, 0.009, -0.007, -0.007, and -0.007, 
respectively. As is seen, the sign of exponent a alternates up to four-loop order and then 
stays negative indicating the stability of the Bose fixed point. Since two neighbouring fixed 
points can not be stable simultaneously we conclude that the fixed point 8 should be unsta- 
ble within the four-loop and higher-order approximations. Consequently, the next, four-loop 
RG approximation will be sufficient to yield the correct structure of the RG flow diagram 
as well as high-precision numerical estimates for the critical exponents. 
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FIGURES 

FIG. 1. Three-dimensional flow diagram of the RG equations ior N = 2 in the three-loop 
approximation. 
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TABLES 

TABLE L Coordinates of the fixed points of the RG equations for N = 2 and = 3 obtained 
within three- loop ([3/1]) and two-loop ([2/1]) approximations. 







1 


2 


3 


4 


5 


6 


7 


8 


N = 2 


Uc 


[3/1] 


0.0 


1.3671 


0.0 


0.9254 


0.0 


0.0 


0.9255 


0.1872 




[2/1]* 


0.0 


1.4863 


0.0 


1.0166 


0.0 


0.0 


1.0166 


0.0332 


Vc 


[3/1] 


0.0 


0.0 


1.8883 


0.7764 


1.6833 


0.9442 


0.3882 


1.4914 




[2/1]* 


0.0 


0.0 


2.1289 


0.8358 


1.8700 


1.0645 


0.4178 


1.8344 


Zc 


[3/1] 


0.0 


0.0 


0.0 


0.0 


1.6833 


2.8325 


1.1646 


1.4914 




[2/1]* 


0.0 


0.0 


0.0 


0.0 


1.8700 


3.1934 


1.2537 


1.8344 


N = 3 


Uc 


[3/1] 


0.0 


1.3310 


0.0 


0.6005 


0.0 


0.0 


0.6005 


0.0780 




[2/1]* 


0.0 


1.4262 


0.0 


0.6307 


0.0 


0.0 


0.6307 


0.0094 


Vc 


[3/1] 


0.0 


0.0 


2.2030 


1.4971 


1.9639 


1.1015 


0.7485 


1.8845 




[2/1]* 


0.0 


0.0 


2.4837 


1.6855 


2.1816 


1.2419 


0.8428 


2.1716 


Zc 


[3/1] 


0.0 


0.0 


0.0 


0.0 


1.9639 


3.3045 


2.2456 


1.8845 




[2/1]* 


0.0 


0.0 


0.0 


0.0 


2.1816 


3.7256 


2.5283 


2.1716 



* Quoted from Ref. § 
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TABLE II. Critical exponents for iV = 2 (upper lines) and N = 3 (lower lines) calculated 



within three-loop approximation. 





i 


o 
z 


o 

O 


4 


r 



D 


i 


Q 

o 


7 


1 
1 






i.oyoy 


i.ouyo 




i.oyyu 


l.OODU 




1 


1.5164 


1.2406 


1.4080 


1.3098 


1.2406 


1.4080 


1.3291 


V 


0.0 


0.0257 


0.0246 


0.0258 


0.0260 


0.0246 


0.0258 


0.0261 




0.0 


0.0238 


0.02 If) 


0.0253 


0.0200 


0.0216 


0.0253 


0.0201 


a 


0.0 


-0.1669 


0.1159 


-0.1258 


0.0094 


0.1159 


-0.1259 


-0.0305 




0.0 


-0.3020 


0.1159 


-0.1390 


0.0094 


0.1159 


-0.1390 


-0.0199 
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